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Abstract 

We examine the possibility that there may exist a logarithmic correction 
to the infrared asymptotic solution with power behavior which has recently 
been found for the gluon and Faddeev-Popov ghost propagators in the Landau 
gauge. We propose a new Ansatz to find a pair of solutions for the gluon and 
ghost form factors by solving the coupled Schwinger-Dyson equation under a 
simple truncation. This Ansatz enables us to derive the infrared and ultraviolet 
asymptotic solutions simultaneously and to understand why the power solution 
and the logarithmic solution is possible only in the infrared and ultraviolet 
limit respectively. Even in the presence of the logarithmic correction, the gluon 
propagator vanishes and the ghost propagator is enhanced in the infrared limit, 
and the gluon-ghost-antighost coupling constant has an infrared fixed point (but 
with a different (3 function). This situation is consistent with Gribov-Zwanziger 
confinement scenario and color confinement criterion of Kugo and Ojima. 
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1 Introduction 



Recent investigations for the infrared (IR) dynamics of Yang-Mills theory based on 
various methods converge to a consistent picture for the gluon and Faddeev-Popov 
ghost propagators: In the framework of Euclidean SU{Nc) Yang-Mills theory with 
the manifestly Lorentz covariant gauge fixing, the transverse gluon propagator in the 
Landau gauge vanishes in the IR limit — > 0, 

hm Dt{p') = 0, (1) 

while the Faddeev-Popov ghost propagator is enhanced in the IR limit: 

hm Gghip') = oo, (2) 

where Dxip^) and Gghip^) are respectively defined by the gluon propagator in the 
Landau gauge and the ghost propagator with unbroken color symmetry as 

Dt^{p) ■■=^^^Dt{p')P^,M. (3a) 
G^^(p) :=5^^G,,(p2), (3b) 

with A, S = 1, ■ ■ ■ , N^ — l and the transverse projection operator, Pj^^ip) := S^^ — ^^^- 

Such behaviors of the gluon and ghost propagators were first derived by Gribov 
long ago as a result of restricting the region of functional integration over the gluon 
field to the interior of the Gribov horizon in order to avoid Gribov copies. Indeed, he 
obtained the result: 

where M is a mass scale called the Gribov mass. Subsequently, Zwanziger has 
obtained the exact analytic result that the restriction of the Gribov horizon enforced 
by a horizon condition yields 

lim Dt{p^) = 0, lim p^Gghip^) = oo. (5) 

p2— >0 p2— +0 

Recent studies of the coupled SD equation for the gluon and ghost 

propagators have reached a conclusion that, in the IR limit — 0, the gluon form 
factor -F(p^) and ghost form factor G{p'^) defined by 

Dt{p') := F(/)/p^ G,.(p2) := G{p')/p\ (6) 

exhibit the power behavior characterized by an IR critical exponent k: 

F{p') = A{py^, G{p^) = B{p')-\ (7) 

where A and B are p- independent constants, and k takes the value between 1/2 and 
1 depending on the approximations adopted to write down the solvable SD equation 
(Gribov's result corresponds to n = 1. However, the aim of this paper is not to obtain 
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the precise value of k.). These IR asymptotic solutions should be compared with the 
ultraviolet (UV) asymptotic solutions obtained by perturbation theory: 

F{p') = Cilogp'r, G{p') = D{\ogpy, (8) 

where C and D are p-independent constants and, at one-loop level, the exponents 
7 and 6 are independent of the number of colors Nc in SU{Nc) Yang-Mills theory, 
i.e., 7 = — ^ and S = This type of solutions with logarithmic dependence on 

momenta can also be derived as the UV asymptotic solution to the SD equation, 
see 1,1. 

It is possible to explain how Gribov's old result is compatible with the recent IR 
solution of SD equation. Gribov's prescription to cut off the functional integral at 
the first Gribov horizon does not alter the SD equations of Faddeev-Popov theory, 
as pointed out by |T^. The cut-off at the first Gribov horizon assures that both the 
gluon and ghost propagators are positive, and that negative and oscillating solutions 
are excluded. Indeed, the above IR solutions for the truncated SD equations are 
positive A,B > 0, 

Moreover, the method of the stochastic quantization in its time- independent for- 
mulation indicates qualitatively the same result |1O|,0 as the SD equation. Recent nu- 



merical simulations on the lattice also support the above results, see e.g., [^,|TB|,|T^,|T3| . 

The above result has an implication to the confinement problem. As an immediate 
consequence of the suppression of the gluon propagator in the IR limit, gluon con- 
finement is easily understood, since the would-be physical transverse gluon becomes 
short range and the massless gluon is absent from physical spectrum. On the other 
hand, infrared enhancement of ghost propagator suggests that the unphysical longi- 
tudinal gluon becomes long range and mediates the confining strong force, providing 
the linear confining potential, i.e., quark confinement. More generally, it is worth re- 
marking that the enhancement of the ghost propagator is consistent with a sufficient 
condition for color confinement (confinement of all the color non-singlet objects) due 
to Kugo and Ojima ||16| (see 



lim G-^p"^) = 1 + u{0) = 0. (9) 

Thus the precise information of the IR limit is quite important to elucidate the non- 
perturbative dynamics of Yang-Mills theory. 

In this paper, we re-examine the truncated Schwinger-Dyson equation for the 
gluon and ghost propagators in Yang-Mills theory in the Landau gauge, keeping these 
developments in mind. A purpose of paper is to point out that the above solution 

(0) is not the most general IR asymptotic solution of the truncated SD equation. 
Another purpose is to understand why the quite different asymptotic solutions, i.e., 
IR power solution and UV logarithmic solution in momentum space, can be obtained 
from the same SD equation. For this end, we study the IR and UV asymptotic 
solution on equal footing simultaneously by proposing a new Ansatz and examine a 
possible logarithmic correction to the IR asymptotic power solution. 

In the studies of SD equation so far , an Ansatz for the solution, (0) or 

(1) , was substituted into the SD equation and the power index and the coefficient are 
determined self-consistently by comparing both sides of SD equation. However, this 
procedure does not guarantee the uniqueness of the solution. Moreover, improving 
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the IR asymptotic solution was attempted by von Smekal, Hauck and Alkofer |Q 
and Atkinson and Bloch 0. It has been shown that the truncated SD equation 
with bare vertex functions is satisfied by the asymptotic expansion (according to the 
notation HI): 



/ TV \ / N \ 



n=l / \ n=l 



(10) 



where k = 0.769479 (p = 1.96964 and /i = l,gi = 0.829602) and a is a free pa- 
rameter which can not be determined from the IR asymptotic expansion alone. The 
IR parameter a is expected to be negative for the IR asymptotic solution to connect 
into the UV asymptotic solution (in a very narrow transition region). The value was 
determined by numerical solution to be negative, a = —10.27685, as expected. See 
for further information. 

In this paper, we adopt another Ansatz for asymptotic solutions which incorpo- 
rates possible logarithmic dependence and can be applied to both IR region — > 
and UV region p"^ —> oo simultaneously: Q 



F{p') = Ae-^'-z^ f 1 + E ^1 , G{p') = Be^^'-/ f 1 + E (12 

\ n=l ^ / V n=l ^ ) 



with a new variable [18], 



p2 

c^ln h-^cr, (13) 



a 

where uj is an appropriate real number and o is the renormalization group invariant 
momentum scale defined by a := /i^exp[— 2 /^^ -^^y] for the renormalization scale /i 

using the /? function: /5((?) := P^^t^- Here the coefficients and dn can be expressed 
in terms of 7 and 5. The IR and UV asymptotic region corresponds to \z\ 00 
where the signature of z depends on that of lo. Note that these parameters are not 
all independent (We can give a certain value for z^, e.g., Zo- = 1 or z^- = 0, since it is 
negligible when | ln^| — > 00). 

An advantage of this Ansatz is that we can search for the possible IR and UV 
solutions simultaneously without discriminating the IR and UV region by specifying 
possible set of power indices and coefficients. In fact, we show that the above Ansatz 
satisfies the truncated SD equation [0 specified below, if the exponents satisfy either 
of the two relationships, 

a = = /3, 7 + 25 = -I,---, (14a) 
a = -2/3 ^0, 7 = -2(5, ■ • • , (14b) 



^In the IR region, we can add the power-series correction part as in (]To|). Moreover, we can also 
treat the region z ~ ^ 1 (j? = a) by changing the last series as 



(11) 
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where the relationships among c„, dn are given later, since they are very complicated 
(see Appendix). It should be remarked again that these relations can be obtained 
without specifying the momentum region in question, by comparing the both sides of 
the respective SD equation. We must choose the appropriate solution as follows. 
The asymptotic UV solution is selected by requiring 

F(p2),G'(p^) ^0 as p2 ^ oo. (15) 

Therefore, the first choice ( |14a|) of parameters corresponds to the UV asymptotic 
solution: pure logarithmic solution without leading power behavior. Indeed, the 
perturbative result 7 = — ^ and S = —-^ satisfies •y + 2S = —1. 

Another advantage is that the IR asymptotic solution allows the logarithmic cor- 
rection to the power solution. It is shown below that the IR solution is identified with 
the second set of parameters ( |14b|) (k = — /5 ^ 0): 

N N 

F{p^) = Ae'""/";^-^^ J2 CnZ-'\ G{p^) = 5e-""/"/ ^ c/^z"", (16) 

n=0 n=0 

where cq = 1 = do and the coefficient c„, dn can now be expressed by 6. Irrespective 
of the value of 6, these solutions for k > behave as 

F(/) ^ (p2)2.(_ iogp2)-25 ^ 0, G{p') ~ {pY^{- logp'Y oo. (17) 

For /t > 1/2, in particular, we have the IR suppressed vanishing gluon propagator 
and IR enhanced ghost propagator, irrespective of the presence or absence of the 
logarithmic correction. Therefore, the existence of logarithmic correction does not 
change the expected IR physics at the IR limit p^ ^ 0. Thus the new type of 
solution, even if it exists, does not destroy the argument in favor of color confinement 
mentioned above. If 5 = and c„ = = (i„ for all > 1, the solution in question 
reduces to the pure power solution known so far. 

The identification of the IR solution is done as follows. First of all, it should 
be remarked that all the solutions of SD equations are not allowed as physically 
meaningful solutions. To clarify this point, we examine the renormalized (running) 
gluon-ghost-antighost coupling constant defined by 

g'{p):=g'F{p')G'ip'), (18) 

where g'^{p) is renormalization group invariant in the Landau gauge, since ZgZ^^'^Z^ = 
Zi = 1. The above Ansatz yields the running coupling g'^{p) = (^2^2^2^(0+2/3)2/1^ ^7+25 

X (^1 + J2n=i CnZ~"'^ (1 + J2n=i ^nZ'"^ ■ The first sct of parameters leads to 

g\p) = g'ABh-' (l + (ci + 2d,)z-' + O {z~')) , (19) 

which vanishes in both UV — > 00 and IR limit — > 0. To be consistent with 
the UV asymptotic freedom, the running coupling must become zero at p = 00 and 
grows as decreases. Therefore, the running coupling (|19]) can be identified with the 
UV solution. However, if we adopt this solution also as the IR solution, the running 
coupling must vanish at p^ = again. This identification leads to the multi-valued (3 
function. Hence, the first set can not correspond to the IR solution. 
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Figure 1: The jS functions l3{g'^) with an infrared fixed point at g^{0) corresponding 
to the IR power solution in the absence (left) and presence (right) of logarithmic 
correction. 

On the other hand, the second choice 

g\p) = g'AB' (l + (ci + 2dr)z-' + O {z-')) (20) 

converges to a finite and non-zero value in both UV ^ oo and IR limit ^ 
irrespective of /5, 5. Hence this is not accepted as the UV solution. Identification with 
the IR solution implies that (7^(0) = g^AB^ is nothing but the IR stable fixed point for 
arbitrary values of k, 5. As p increases towards the UV region from the IR fixed value 
at p = 0, therefore, the running coupling must bend below to smoothly connect into 
the UV logarithmic solution. This effect can be seen by taking into account the power 
series in . For this purpose, the parameters ci and di must satisfy {ci + 2di) / z < 0, 
provided that cirfi ^ 0. Incidentally, the IR diverging running coupling constant does 
not appear from the solutions subject to the above relationship of parameters. 

The UV behavior is characterized by the j3 function, j3{g'^) = — ^^^j/^ g'^+o{g^) with 
P'{0) = 0. In the neighborhood of IR fixed point g'^{0), the beta function Pi^g"^) := 
jJ'^^^-J^ behaves as follows (See Fig. |l|). The power solution with power correction 
(Pleads to P{g') = 2pg{0ngyg\0) - 1) + o{{gyg\0) - 1)^), and (^'{g'm = 
2p > 0, while the power solution with logarithmic correction ([TBD yields l3{g'^) = 
-^^(^/V^/'(0) - iy + o{{gyg\0) - 1)3), and p'{g\0)) = 0. The difference can not 
be seen in the g'^{p) vs. \np plot, since = at p = in both cases. However, 

in the g^{p) vs. p plot, the logarithmic correction implies = —oo at p = 0, in 

contrast to the power solution ^^^^^ = at p = 0. The data of numerical simulation 
available seem to be still insufficient to discriminate between two cases. In this paper, 
therefore, we examine the IR behavior in the analytical way. 

2 Finding asymptotic solutions of the coupled SD 
equation 

For simplicity of analysis, we restrict our consideration to a version of the truncated 
SD equation adopted by Atkinson and Bloch using the bare vertex functions, 
since qualitatively the same IR behavior was obtained also from the solutions of the 
truncated SD equations using the full vertex function improved so as to be consistent 
with the Slavnov- Taylor identity. 

The truncation is done as follows. We include all diagrams in the ghost equation, 
while we neglect contributions from the two-loop diagrams in the gluon equation. The 
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Figure 2: Diagrammatic representation of the truncated Schwinger-Dyson equations 
for gluon and ghost propagators. 



gluon equation was contracted with the Brown- Pennington projection operator |20 
Rfiu{p) '■= — 4^p^, which removes the tadpole term. Therefore, the quadratic UV 
divergence can be removed by this projection. Finally, the non-trivial contribution 
in the gluon equation comes from only two diagrams, i.e., gluon loop and ghost loop, 
apart from the tree gluon propagator. We include only the bare triple gluon vertex 
and the bare gluon-ghost-antighost vertex. See Fig. ^. 

Moreover, in order to avoid angular integration, we adopt the approximation 
called the ?/-max approximation, — g)^) = F(max{p^, g^}), G{{p — qY) = 

G(max{p^, g^}). Here the angle 6 comes from the inner product p ■ q = cos6 
between the external momenta p and the internal momenta q which is to be integrated 
out. Thus we can write explicitly down the coupled SD equation each of which is the 
one-dimensional integral equation. 

Introducing the renormalization constants for the gluon field Z^, the ghost field 
Z3, the triple gluon vertex Zi and the gluon-ghost-antighost vertex Zi, we write down 
the SD equation for the renormalized gluon and ghost form factors. By introducing 



A : = 



167r2 ' 



X 



and y := g^, the ghost SD equation reads 



G-\x) 



3iV, 
4 



\Zi 



Fix) 



dyyG{y) + 



dy 



F{y)G{y) 



(21) 



while the gluon SD equation reads 



F-\x) 



i 3 



G{x) 



3 



lF{x) 



dyy'Giy) + ^ 

2x^ 
17F{x) 



dyyG{y) + \f'^' ^G^y) 
2 Jx y 



Jx y 8 Jx 



dyyF{y) - 



9F(x) 



5X 



dyF{y) 



^This approximation was also called Landau- Abrikosov-Khalatnikov (LAK) approximation or 
Higashijima-Miransky (HM) approximation. In the bare vertex approximation, the exact angular 
integration was performed for the pure power solution |^ without relying on this approximation, 
which yields a largest value of 1. This value agrees with the Gribov prediction. Recently, 

it has been shown [|[|| that the value of n is decreased between 0.5 and 0.6 due to the renormalization 
effect of the triple gluon interaction vertex which is necessary to reproduce the correct coefficient of 
/3 function in the UV limit. 
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The renormalization constants Z3 and Z3 are eliminated by subtracting the equation 
at X = a. The ghost equation reads 



3N, 



^^""^ dyyG{y) - [x ^ a) + f '^F{y)G{y) 

Jx y 



x^ Jo 



(23) 



while the gluon equation reads 



F~Hx) =F-Ha) 



Nr. 



'-\Zi 



a;3 Jo 



dyy'G{y) + 



3G(x) 
2x2 



dyyG{y) 



[X cr 



3 

7x 



dyy'F{y)- 



2x3 



2x2 



dyyF{y) 



9F(x) 



5X JO 



+ ^ ^F\y)-{x^a)-l ^F\y) 



X y^ 



dy 



X y 



dyF{y) 
(24) 



First, we rewrite the coupled SD equation in terms of the new variable z 
+ uj\n.- and C := z^- + uj In-. The ghost equation reads 



G-\z) - G-\z^) 

and the gluon equation reads 

F-\z)-F-\z^) 



3iV, 



F{z)e-^'l^ r ^e^i/^GiO 

J—ujoo CO 



^^^a;z^z^)+ ^F(C)G(C) 

UJ 



dC 



(25) 



G{z)e 



^e^^l-G{Q + hi{z)e-^^'- 



[X a;z ^ z„) + 



LJOO UO 

dC 



dC 



woo id 



e^^'-G{0 



UO 



G\0 
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2 J— woo 2 J—UJOO CO 

r ^e^/^FiO + lei- r ^e-</-F\() 

J— woo CO 8 J z CO 



-{x~^a-z^z„)-7 -F\C) 



UO 



(26) 



Next, we substitute the Ansatz ([T^) into the above equation and perform the 



integration over C, by making use of the integration formula |21 

^ J v-l)-''«-^'+'^r[l + 6, -az] {a ^ 0) 
^ ^ ^'l + hr^z^+' (a = 0) 



(27) 



■^In order to use the asymptotic solution in the integrand, we perform the decomposition: 

dyfiy) = - 1^2 dyfijj) + dyf{y) for the IR hmit x 0) and dyf{y) = J^^ dyf{y) + 

/o dyf{y) for the UV hmit x — > A^(^ 00). In our approach, the term dyf{y) can not be 
determined, but it is a p-independent constant and does not affect the exponent. 
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where r[c, x] is the incomplete gamma function with the asymptotic expansion for 
large \x\: r[c,x] = x^-^e"^ [l + E^=i ^''"^^^''"n^"'^''"''^ ] ■ In particular, for large z and 
a 7^ 0, we can use the formula, 



1 + E 



(28) 



The ghost equation is integrated out, when a ^ —P and /3 7^ —2: 



4 

g(a+/3)2:A 

a + P 



2 + /3 ' 



a;5 



+ ci^^-^ + o(^^-')) <^ / - ^7^^ + diz'-^ + o{z 

L 2 + p 



(29) 

where we have required /3 > — 2 for eliminating the terms coming from the lower 
bound —ijjoo to avoid the IR singularity. By comparing both sides of this equation, 
we find that the following relations must be satisfied. 



a + 2/5 = 0, 
7 + 25 = 0, 



1 ^ ^AZiAS^ 



1 1 
+ 



Ci + — Ci + (ii — 



(30a) 
(30b) 

(30c) 
(30d) 



2 + /3 « + 

Similarly, the gluon equation reads, if o; 7^ 1/2, 0, —1, —2, —3 and /9 7^ 0, —2, —3: 

^-lg-a./a;(^-T _ ^^^-7-1 + o(C'^-^)) - {{z ^ Z,]) 



3 

3 e^^^l^ 



3 + /?' 



+ di/-^ + o(C'^"')) <! 



cj5 
3 + /?' 



J-i 



(/ + d,z'-' + o(C^^-2)) |/ - ^z'-' + diz'-' + o(C^^-2)| 



.2a 



- ^^^25-1 ^ 2rf,;,2a-i + o(C^^-2)| - {[z ^ ;.,]) 



/5 

--^ ^2az/u 



23 + a 



{z^ + c,z^-' + oiC^-')) [z-^ - -^z^-' + cz-^-' + o(C^^-^)} 



8 l + a 

7 e2o:./. 



(z^ + CiZ^-^ + o(C"^-2))|^^ 
(^^ + Ci^^-^ + o(C'^-2))|z^ 



2 + a 



3 + q; 



-^ + o(C^n} 



2 2 + a 



+ 7- 



^2a- 1 

g2aC/^ 



.27 



2UJ' 



^^27-1 



2a 



2a -I 

27 _ ^^27-1 + 2ciZ^^-^ + 0(C'^-')i - ([^ - ;^.]) 



(31) 



where we have required a > —2 and (3 > —2 for ehminating the terms coming 
from the lower bound — cjoo. In the right hand side of (pT]), the former half is the 
contribution from ghost loop and the latter one comes from gluon loop. If a > 0, 
then j3 = —a/2 < and hence the ghost loop is dominant over the gluon loop. This 
phenomenon was called the ghost dominance in the pure power solution. We confirm 
that the ghost dominance also holds even in presence of logarithmic correction (under 
the truncation adopted). 

The parameters must satisfy the following relations. 



a + 2/3 
1 + 26 



0, 
0, 

3 

ci + di^XZiAB 

= ^^XZ.AB' 
3 




(32a) 
(32b) 

(32c) 



(3 + /5)2 2(2 + /?)2 4/?2 



(32d) 



The relation ( |30a|) is equal to ( P2a| ), and the relation (|30b|) is equal to (P2b|) . From 
the relation (pOc|) and (|32c|), we obtain 



-1 



XZiAB^ 



+ 



3 + P 2 + 13 

Therefore, [3 satisfies the algebraic equation, 19/5^ + 77/3 + 48 
are /? = Tl±^^ ^ -0.769479, -3.28315, irrespective of N, 




(33) 



0, whose solutions 
By adopting (3 = 

— 0.769479 (> —2), the IR fixed point of strong fine structure constant is given by 

a,(0) := AnXZiAB'^ = {g'^ / A'k)ZiAB'^ = (3/iV,)11.4702. (34) 
The running coupling in the neighborhood of = is given by 

c^sip) = {g^/A'K)ZiF{p')G\p) = (3/iVe) 11.4702 [l + (ci + 2di)z~^ + O [z'^)] . (35) 

The relationships among c„, dn {n > 2) are given in Appendix where trick ending 
results 7 = = 5 and c„ = = (i„ are obtained from the consistency of the coupled 
SD equation. 

Another solution is obtained as follows. Both the gluon and ghost loops contribute 
equally to the leading order, only when a = = (3. From the matching of the last 
term in the respective contribution, we obtain 

= p, and -7 = 2(5 + 1 = 27 + l,i.e.,7 = -1/3 = 5. (36) 



a 



This is the same result as the UV case of Atkinson and Bloch [Q. If a = 0, then 
(3 = 0, and vice versa. This solution implies the UV asymptotic freedom, although 
it leads to a different value Po = ^NcZi for the coefficient of the /3 function from the 
perturbative result Pq = ^-NcZi where Zi = 1. This disagreement was resolved by 
taking more involved renormalization prescription in The correction terms with 

coefficients c^, dn can be obtained in a self- consistent way. The more details will be 
given in a subsequent paper. 
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3 Conclusion and discussion 



We have proposed a new Ansatz for the asymptotic solution of the coupled SD equa- 
tion for the gluon and ghost form factors in Yang-Mills theory. This Ansatz can be 
apphed to find the UV and IR asymptotic solution simultaneously. In fact, we have 
reproduced the UV asymptotic solution with logarithmic behavior and IR asymptotic 
solution with power behavior. Moreover, we have pointed out a possibility that an IR 
asymptotic solution can have a logarithmic correction to the power behavior found 
so far, with keeping the IR fixed point. However, it is proved in Appendix that the 
logarithmic correction in the IR asympotic solution is excluded for a version of the 
truncated SD equation treated in this paper. In this case, therefore, ( |14b|) is replaced 
by 



a 



-2/5^0, 7 = = 5, c„ = = 4 



(37) 



and only the power correction ( p^O]) is allowed in the IR region, which was presumed 
in the previous papers . This result is in sharp contrast with the UV asymptotic 
solution which allows power corrections to the logarithmic behavior . 
The issues to be clarified in subsequent works are as follows. 

1) We are to study the renormalization group properties associated with the IR 
asymptotic solution, e.g., the /3 function and the anomalous dimension near the IR 
fixed point. They might be helpful to exclude the logarithmic correction in the IR 
asymptotic solution. The result will enable us to select the physical solution with 
renormalization group invariance, since all of the solutions of SD equation are not 
necessarily the physical solution. 

2) In this paper we have adopted the bare vertex function. Therefore, it will be 
important to see whether our result remains true even after the vertex function is 
improved so as to satisfy the Slavnov- Taylor identity. 

3) It will be interesting to study the IR asymptotic solution in other covariant 
gauges, e.g., the Maximal Abelian (MA) gauge pS), in addition to the non-covariant 
gauge, e.g.. Coulomb gauge p3|,^]. We wish to know how the behavior of IR solution 
in the MA gauge is connected to the dual superconductor picture 
confinement 



26[| for quark 



27 



4) We have not yet worked out the effect of the dimensionality of space-time on 
IR solution. For example, we can ask whether the presence or absence of the power 
or logarithmic correction is specific to the four-dimensional case or not. 

5) It is not yet clear how the UV solution with gluon condensation [^,|2S|,^ is 
connected to the IR solution in a quite narrow transition region. Does the mass scale, 
e.g., the Gribov mass in IR region have something to do with the gluon condensation? 
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A Proof of absence of logarithmic corrections in 
the IR limit 

Substituting the Ansatz ([T2|) into the SD equation and performing the integration 
exphcitly, we obtain the recursion relations by equating the coefficients of in 
both sides: 



n,m>0:n+m=N 



SAT, 



N 

E 



1=1 ^K'^ ^ yj n,m>0:n+m=N-i 

E dmCn\{{m + n-^ - 5 + (38a) 



1=0 

e-i 



n,m>0:n+m=Af 

AT, ^ 



-u 



E 

£=1 



^ n,m>0:n+m=Ar-£ 

dm.d.^ 

3 



+ 



UJ 



(3 + /?)^+i 2(2 + /?)^+! 

1 

2W) 



J n,m>0:n+m=N-£ i=0 



e+1 



1 

E dmdnY[{m + n-26 + 

n,m>0:n+m=N-£ i=0 ^ 



(38b) 



where Gn and Fn are defined by ^-^(2) = A-^e-'^''^ z~^' 'En=o Fnz'^ and G'^z) 
B~^e~^^^'^z~^ Y.n=oGnz~'^ and exphcit calculations show 



Gi — — di, G2 — —d2 + di, G3 
G4 = — 6/4 + 2did2, — ?>d\d2 + d\,- ■ 



Cl, 



-C2 + q, 



-4 + 2did2 - rfi^, 
-C3 + 2CiC2 — c'^, 



-F4 = - C4 + 2CiC3 - 3c^C2 + cf, 



(39) 



It is easy to find that a + 2/5 = and 7 + 25 = 0. 

In the following, we shall prove 7 = = 5 and Cjv = = rfjy for A^ > 1 by 
mathematical induction. First, we show that 7 = = 5. For A^ = 1, the above 
recursion relations lead to 



-di - {di + Cl) = - 



3iV, 

4 
A^. 



-V 



+ 



(2 + /3)2 (a + /3)2 
1 3 1 



o;(-5), 
1 



"""'^ Y''l(3 + /3)2 2(2 + /3)2 ' (2/3)2 



+ 



. Eq.i 



and 



(40a) 

cj(-5), (40b) 
agree with (^Od|) and 



where we have used 7 = —25 to derive 
( P2cl| ), respectively. Substituting a = —2(3 and [3 = —0.769479 into the RHS leads 
to different coefficients of AcZ/co'5, although the LHS are the same, i.e., Ci + 2di = 
— 1.762 A"cZ/C(j5 and 0.12247Aci^co'5. Therefore, we conclude 5 = and hence 7 = 0. 
Thus we obtain ci + 2di = or ci = —2di. 



11 



For = 2, we use 7 = = 5 to write the recursion relations, 
-d, + d,- {d, + d,c, + c) = - + (41a) 

2 /, ,2 , N r 1 3 1 11, 

-c, + c, - (d, + rf, + d,) = - - 2 + j-^^- (41b) 

For Ci = — 2(ii, both LHS agree with — C2 — 2^2 + Bc?^. However, the coefficients of 
Nf.i'ujdi in RHS are not equal to each other. Hence, we obtain di = and Ci = 0. 
Thus, we obtain C2 + 2^2 = 0. 

For N >3, suppose that = = rf^ for m = 1, 2, ■ ■ ■ ,N — 2 and ciy_i + 2dN-i = 
(Note that cq = 1 = do). Then, the non- vanishing contributions in RHS of (|38a]) or 
( P8b|) stem from the terms containing the non- zero c^^i or d^^i. Hence it is sufficient 
to consider the i = 1 part in RHS. Thus, the recursion relations reduce to 

Gn - {dN + Civ) = - -^^^1 /2+^2 ^^-1 - -(—]^{dN-i + c^-i) j(A^ - 1) 



/o , m2 + , r - (42a) 
^ - 2c?jv = - TTT^m? - + TTT^ - l)c?iv-i^. (42b) 



(2 + /3)2 (a + /3) 

3iVc f 1 1 

4 ''l(2 + /3)2 + (a + /?)2 

3''l(3 + /3)2 2 (2 + /3)2 + (2/3) 

Under the assumption, we have the simple expression, = —d^ and F/v = —cjq. 
The LHS of both equations reduce to —(cat + 2(iiv). By substituting the value of 
f3 = —K into the RHS, we find that the respective RHS gives the different coefficient 
for Ncui^N — l)djq_iuj. This implies that d^-i = and hence CAr_i = 0. Thus 
we obtain cn + ^d^ = 0. Repeating this argument, we obtain = = dj^ and 
cn+1 + 2dN+i = 0. Thus we can conclude cat = = c^at for all > 1. 
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